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Abstract. We review recent results on the properties of materials with corre- 
lated electrons obtained within the LDA-I-DMFT approach, a combination of a 
conventional band structure approach based on the local density approximation 
(LDA) and the dynamical mean-field theory (DMFT). The application to four 
outstanding problems in this field is discussed: (i) we compute the full valence 
band structure of the charge-transfer insulator NiO by explicitly including the 
p-d hybridization, (ii) we explain the origin for the simultaneously occuring metal- 
insulator transition and collapse of the magnetic moment in MnO and Fe2 03, (iii) 
we describe a novel GGA-I-DMFT scheme in terms of plane- wave pseudopotentials 
which allows us to compute the orbital order and cooperative Jahn- Teller distor- 
tion in KCuFa and LaMnOa, and (iv) we provide a general explanation for the 
appearance of kinks in the effective dispersion of correlated electrons in systems 
with a pronounced three-peak spectral function without having to resort to the 
coupling of electrons to bosonic excitations. These results provide a considerable 
progress in the fully microscopic investigations of correlated electron materials. 

1 Introduction 

The theoretical understanding of materials with strong effective interactions between the elec- 
trons is one of the most challenging areas of current research in condensed matter physics [I] . 
Strong electronic interactions are encountered in materials with open d and / shells, such as the 
transition metals Ti, V, Fe and their oxides or rare-earth metals such as Ce, where electrons 
occupy narrow orbitals. This spatial confinement enhances the effect of the Coulomb interac- 
tion between the electrons, making them "strongly correlated" . The interplay between the spin, 
charge and orbital degrees of freedom of the correlated d and / electrons and the lattice degrees 
of freedom leads to a multitude of unusual ordering phenomena at low temperatures. Conse- 
quently, strongly correlated electron systems are often exceedingly sensitive to small changes in 
the temperature, pressure, magnetic field, doping, and other control parameters. This results, 
for example, in large changes of the resistivity across metal-insulator transitions or upon the 
application of a magnetic field of the volume across phase transitions, and of the effective elec- 
tronic masses; electronic correlations are also essential for an understanding of high-temperature 
superconductivity. Correlated electron materials often reveal rich phase diagrams originating 
from the interplay between electronic and lattice degrees of freedom |ll2l3j . These compounds 
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are particularly interesting in view of possible technological applications. Namely, their prop- 
erties can be employed to construct materials with new functionalities for future technological 
applications. 

The electronic properties of many materials can be computed using density functional theory 
methods, e.g., in the local density approximation (LDA) [4l, the generalized gradient approx- 
imation (GGA) |5|6| . or using the so-called LDA-fU method 7 8 9. Applications of these ap- 
proaches often describe the physical properties of simple elements and intermetallic compounds, 
and of some insulators, quite accurately. Moreover, they permit to make correct qualitative pre- 
dictions of the magnetic, orbital, and crystal structures of solids where the equilibrium (thermo- 
dynamic) structures are determined by simultaneous optimization of the electron and lattice 
systems [10 11 12 13^. However, these methods usually fail to describe the correct electronic 
and structural properties of electronically correlated paramagnetic materials. Hence the com- 
putation of electronic, magnetic, and structural properties of strongly correlated paramagnetic 
materials remains a great challenge. 

Here the recent combination of conventional band structure theory [4 and dynam- 
ical mean- field theory [14115116117118) . the so-called LDA-I-DMFT computational scheme 
[19120121122123124] has led to a powerful new tool for the investigation of strongly correlated 
compounds both in their paramagnetic and magnetically ordered states. In particular, it has 
already provided important insights into the spectral and magnetic properties of correlated 
electron materials [25l:^6i:j7l:j8l^9l30l3ll32l33l34l35l36l37l38l39l40l41l42l43l44] especially in the 
vicinity of a Mott metal-insulator transition as encountered in transition metal oxides. 

In this paper we review several topics related to the physics of correlated electron materi- 
als where significant progress was made during the last three to four years. In particular we 
describe the application of the LDA+DMFT approach to charge-transfer materials where — 
in contrast to the early transition metal oxides — additional complexities from the presence of 
p bands arise. We apply the scheme to study the single-particle spectrum of the prototypical 
charge-transfer insulator NiO. By explicitly including the 0-p orbitals and their hybridization 
with Ni-d orbitals we obtain a unified description of the full spectrum (Sec. [5]). A second topic 
is the investigation of the infiuence of pressure on correlated electronic materials and of the 
collapse of the magnetic moment. Indeed, in numerous transition-metal oxides the pressure- 
induced metal- insulator transition is accompanied by a collapse of the magnetic moments, i.e., 
a transition from a high spin state to a state with low spin. Such a transition is observed for 
example in MnO and Fe203. Our results obtained with LDA-I-DMFT for MnO are described in 
Sec. [31 A breakthrough in the application of the LDA-I-DMFT scheme concerning the compu- 
tation of correlation induced atomic displacements and structural transformations is presented 
in Sec.m Namely, by formulating a LDA/GGA+DMFT scheme implemented with plane-wave 
pseudopotentials it has been possible to describe and understand the structural transformation 
of paramagnetic solids due to electronic correlation effects. Finally, we present a general expla- 
nation for the appearance of kinks, i.e., sudden changes in the slope of the effective electronic 
dispersion of the correlated electrons. By a detailed analytic investigation of the Green function 
and the self-energy for correlated systems with a pronounced three-peak spectral function we 
are able to explain the occurrence of kinks without having to resort to any explicit coupling of 
electrons to bosonic excitations (Sec. [5]). 



2 Charge-transfer compounds 

In the so-called charge-transfer (CT) materials the ligand orbitals play an important role since 
they determine the low-energy physics of these systems. The concept of a CT insulator was 
introduced in the mid 1980's by Zaanen, Sawatzky and Allen (ZSA) in their classification of 
transition- metal oxides (TMOs) and related compounds PS^. In the early TMOs the ligand 
p band is located well below the transition metal d band and thus plays a minor role in the 
low energy dynamics. Such a case, called Mott-Hubbard system in the ZSA scheme, is well 
described by a multi-band Hubbard model. On the other hand, the late TMOs are of the 
CT type where the p band is located between the interaction split d band. A more general 
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Hamiltonian where the p states are explicitly included is then needed, which can be viewed as 
a combination of multi-band Hubbard and Anderson lattice models. The results of calculations 
using the static mean-field approximation such as LDA+U I7I8I9I suggested that the treatment 
of the dynamical correlations and the ligand-transition- metal hybridization might be important 
to properly describe the excitations in the CT materials. 

NiO is a well-studied prototype of a CT insulator. It is a type II anti-ferromagnet (T^r 
= 523 K) with a magnetic moment of almost 2^b and a large gap surviving well above T/v. 
Its photoemission spectrum |46l47l48l49j . which is our main concern, has several characteristic 
features: (i) a broad high-energy peak with little dispersion dominated by Ni d character, (ii) 
a flat low-energy band with substantial amount of Ni d character, (iii) dispersive O p bands at 
intermediate energies, and (iv) a charge gap of 3-4 eV. All of these feature are essentially the 
same in both the antiferro- and the paramagnetic phase [50] ■ Although further details can be 
resolved in the spectra of NiO we will focus on these four features here. 

The LDA-I-DMFT calculations proceeds in two independent steps: the construction of the 
effective Hamiltonian from converged LDA calculation, and the subsequent solution of the 
corresponding DMFT equations. We use the Wannier function basis ^51l52j for an eight-band 
p-d Hamiltonian 

+ E Urp'nt^nf^,,. (1) 

Here dfcao- and pkja are Fourier transforms of diaa- and Pija, which annihilate the doi p electron 
with orbital and spin indices aa or in the ith unit cell, and nf^^ is the corresponding 
occupation number operator. The elements of the U^^ matrix are parameterized by U and 
J. The constrained LDA calculation yields J7 = 8 eV and J — 1 eV [7 . To account for the 
Coulomb energy already present in LDA we subtract a double counting correction from the 
dd-diagonal elements of the LDA Hamiltonian: 

hi% = hiUik) ~ (iVorb - l)Un^T,K5^p (2) 

where tilda is the average LDA occupation per orbital and A"orb = 10 is the total number of 
orbitals within the shell. With the kinetic and interaction parts of the Hamiltonian specified 
we solve the DMFT equations on the imaginary frequency axis, a key part of which is the 
auxiliary impurity problem treated with the quantum Monte-Carlo (QMC) method |53j. The 
results reported here were calculated at T = 1160 K. To obtain the single-particle spectral 
functions an analytic continuation to real frequencies is performed using the maximum entropy 
method [5^ . 

We start the discussion of our numerical results with the single-particle spectra of stoi- 
chiometric NiO. Using the full p-d Hamiltonian we are able to cover the entire valence and 
conduction bands spectra. In the left panel of Fig. [T] we show the local (fc-integrated) spectral 
densities corresponding to Ni d Cg and t2g, as well as O p electrons. In the inset we compare 
the total Ni d spectral density to the experimental data of Sawatzky and Allen ^46j . The theo- 
retical spectrum very well reproduces the experimental one, including the size of the gap, the 
d character of the conduction band, the broad d peak at —9 eV, the position of the p band, 
and the strong d contribution at the top of the valence band. The dominant feature of the 
valence spectrum is a distribution of spectral weight between the broad peak at high energy 
and sharp peaks at low energy. The spectral weight distribution is a combined effect of the p-d 
hybridization and the electronic correlation in the Ni d shell as was explained by Fujimori et 
al. [SS^ in terms of eigenstates of a NiOe cluster, and Sawatzky and Allen PBI. Emitting a d 
electron from the state the system can reach either the final state (high-energy peak) or 
a d^L_ final state (sharp low-energy peak), with a ligand hole, due to the p-d electron transfer. 

Virtually all electronic structure computational methods have been applied to NiO, often 
claiming success in describing its photoemission spectrum. However, these claims are often 
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Fig. 1. (color online) Left panel: the orbitally resolved single-particle spectrum of NiO with contribu- 
tions from Ni d-Cg (black), Ni d-t2g (red), and O p (blue shaded) states. The inset shows comparison 
of the total Ni d spectral density (with Gaussian broadening of 0.6 eV) (black) to the XPS spectra 
(brown shaded) of Ref. [46] . Right panel: the orbitally resolved single-particle spectrum of hole doped 
NiO (uh = 0.6). The total Ni d spectrum (black) is compared to the experimental photoemission and 
inverse-photoemission spectra of Lio.4Nio.6O [56^ (brown shaded). 



based on picking only a particular aspect of the problem. Putting aside the fact that static 
methods such as LDA-I-U [7I8I9| and generalized density functionals |5I6I as well as GW |57l58j 
provide gapped spectra only in the AFM phase, the main difficulty appears to be capturing 
the distribution of spectral weight between the high-energy dJ peak of Ni d character and the 
low-energy d^L band with mixed Ni d - O p composition. For example, the high-energy Ni d 
peak is completely missing in the calculations with the recently introduced density |5|6j and 
density matrix [SHj functionals and also in GW calculations [57158160161] for reasons explained 
in Ref. [HI]- Also LDA-I-U calculations with small U miss the high-energy peak |63 . On the 
other hand, for larger U values the d character is missing in the low-energy peak |7l61j . The early 
LDA-I-DMFT calculations of Ren et al. [BJ also do not yield qualitatively correct valence band 
since the O p states were not included in the effective Hamiltonian. On the other hand many- 
body techniques such as cluster exact diagonalization [S5] , three-body scattering technique [5S] 
or DMFT ^41.42.66.67.68] , which are capable of generating new poles in the single-particle self- 
energy, describe the distribution of the valence spectral weight in agreement with experiment. 

To further investigate the electronic structure of NiO we have calculated the fe-resolved 
spectra, which determine the correlated bandstructure. In Fig. [2] we compare the theoretical 
bands, represented by the fc-dependent spectral density A(k, w), along the X—F and F—K lines 
in the Brillouin zone with ARPES data of Refs. [48149] . Both theory and experiment exhibit 
two relatively flat bands at —2 and —4 eV followed by several dispersive bands in the —4 to 
—8 eV range and a broad incoherent peak around —10 eV. Overall we find a good agreement. 
The deviations around the F point in the right panel of Fig. [5] are due to the inaccuracy in 
the location of the F point in the off- normal-emission experiment |3H]- The crosses near the 
F point mark a weak band which was interpreted as a consequence of AFM order 149j and 
is therefore not expected to be found in the paramagnetic phase investigated here. Additional 
weak structures in the uppermost valence band (not shown here) were identified in the spectra of 
Ref. [31] . To assess their importance the reader is referred to the original experimental spectra. 

The effect of doping holes to a CT insulator is of prime importance, particularly because the 
high-temperature superconductivity is widely believed to arise from this effect. Besides that, 
the p character of the doped holes is a defining feature of CT materials. In NiO hole doping can 
be achieved by Li substitution. Lia;Nii_a;0 was studied in the doping range x = 0.02-0.4 [55] . 
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Fig. 2. The fc-resolved total spectral function A{k,tjj) along the K ~ F (left panel) and F — X (right 
panel) lines in the Brillouin zone depicted as a contour plot. The symbols represent the experimental 
bands of Shen et al. [ffl]. The theoretical gap edge was aligned with the experimental one. 

Table 1. Orbital occupancies and the local moment on the Ni site for different hole dopings. 
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Starting from the undoped system the replacement of x Ni^+ ions by Li^+ ions introduces 
on average = x/{l — x) holes per Ni site. In the following, we neglect the substitutional 
disorder and study the Hamiltonian of NiO at various non-stoichiometric fillings. As crude as 
this approximation may be we believe that the essential physics of p-d spectral weight transfer 
is captured correctly. In Table [T] we compare the orbital occupancies and local spin moment 
for different hole dopings. While a substantial d spectral density is observed at the top of NiO 
valence band the doped holes reside almost exclusively in oxygen p orbitals. This apparent 
contradiction to the experience from non-interacting electron systems reflects the importance 
of many-body effects. 

Yet another demonstration of electronic correlations is a substantial reconstruction of the 
spectral function, in particular the Cg one, shown in the left panel of Fig. [1] Most notably the 
Mott gap is filled, while the Hubbard subbands remain as distinct features with the spectral 
weight reduced as compared to pure NiO. This is also observed in experiment, as shown in the 
inset. To understand the spectrum it is useful to think about the doping in two steps. First, 
the holes are doped into a system with no p-d hybridization. For the present heavy doping 
this simply amounts to shifting the chemical potential deep into the non-interacting p band. 
Second, the p-d hybridization is turned on and the system is essentially described by the periodic 
Anderson model with a well known spectrum similar to that for the Cg orbital |69) . 

While the d occupancy is barely affected by the hole doping, the behavior of the local spin 
moment changes. The instantaneous moment Ms — \/\>n^ is reduced due to the suppression 
of the dominant d^, iriz — 2 contributions in favor of d^, d^, and d^ contributions with = 
and TTiz = 1. Even more pronounced is the effect of increased screening of the local moment, 
reflected in faster decay of imaginary time spin-spin correlations (see Fig. [3]), leading to a 
significant departure of the screened moment AfcS — \/Tx\oc from its instantaneous value. 

To summarize, we have presented LDA-I-DMFT calculations on NiO and its hole doped rel- 
atives using a Hamiltonian spanning Ni 3d and O 2p bands. It appears that the main limitation 
of the truly first principles character of LDA-I-DMFT studies is the need for a double-counting 
correction, the universal form of which is not known, and perhaps does not exist. We also 
mention the density-density form of the electron-electron interaction used here, a rather techni- 
cal approximation, which can be avoided with the most recent CT-QMC implementations |71| 
or exact diagonalization [7D]. Despite these approximations we were able to obtain the main 
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characteristics of the photoemission and inverse-photoemission spectra of both pure and hole- 
doped NiO. The computational methods claiming success for NiO and CT compounds in general 
should be able to capture substantial d spectral weight at the top of the valence band together 
with the incoherent lower Hubbard band, and, perhaps most importantly, should be able to 
describe the peculiar behavior upon hole doping. 



3 Metal-insulator and spin transitions 

A very interesting type of metal-insulator transition (MIT) is the pressure driven transition 
accompanied by a change of the local spin state (high spin (HS) to low spin (LS) transition) seen 
in MnO [72], BiFeOa [73] or Fe2 03. An understanding of the pressure-driven HS-LS transition 
and its relationship to the MIT and to structural and/or volume changes is relevant to a broader 
class of oxides [74 , often with geophysical implications [7S] . 

We have studied two isoelectronic materials which exhibit the above behavior under pressure, 
hematite (a-Fe203) [31| and manganosite (MnO) [l^. At ambient conditions, Fe203 is an 
antiferromagnetic (AFM) insulator (T^r = 956 K) with the corundum structure [7^ while MnO 
is antiferromagnetic only below 118 K and possesses rock-salt structure. In both compounds 
the transition metal (TM) ions have formal valence state and are octahedrally coordinated 
with oxygen, which gives rise to the characteristic &g-t2g splitting. Photoemission spectroscopy 
(PES) at zero pressure |77I78I79] classified Fe2 03 as a charge-transfer insulator with the charge 
gap of 2.0-2.7 eV inferred from the electrical conductivity data [80 81^. Under pressure, a first- 
order phase transition is observed at approximately 50 GPa (82% of the equilibrium volume) 
with the specific volume decreasing by almost 10% and the crystal symmetry being reduced 
(to the Rh203-II structure) |82I83I84] . The high-pressure phase is characterized by a metallic 
conductivity and the absence of both magnetic long-range order and the HS local moment [82] . 
Badro et al. showed that the structural transition actually precedes the electronic transition, 
which is, nevertheless, accompanied by a sizable reduction of the bond lengths [55] . 

MnO shows a similar behavior. The shock data '5S1, and then Raman and optical stud- 
ies [5Ti8 8i . had identified a transformation in MnO in the vicinity of 90-105 GPa. Transport 89 , 
magnetic [89], structural and spectroscopic |90l91j . and refiectivity ^ data all point to a 
first-order, insulator-metal Mott transition near 100 GPa with volume {v = V/Vq, Vq = vol- 
ume at ambient pressure) collapse w=0.68— >-0.63, and moment collapse (from ^5^b to l/is or 
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Fig. 4. Left panel: the orbital-resolved occupancies of the 3d shell of Fe in Fe203 (red) and Mn in 
MnO (black). The thick lines show the average occupancy Ud per orbital. Right panel: the local spin 
moments in Fe203 (red) and MnO (black) measured as instantaneous moment (diamonds) and screened 
moment (circles). The symbols mark the actual numerical results, the lines are guides to the eye. Vo is 
the volume at ambient pressure. Adapted from [43| . 



less [90l91j ). The structural data indicate a Bl— >B8 change just before the Mott transition, 
which thus occurs within the B8 (NiAs) phase rather than the Bl (NaCl) phase. Since the local 
environment of the Mn ion remains the same, this structural change is not expected to have 
much effect on the Mott transition. 

Our calculations proceed along the lines described in the preceding section: i) construction 
of the multi-band Hubbard Haniiltonian spanning the TM d and ligand p bands with local 
Coulomb interaction within the TM d shell in the density-density (Ising) approximation, ii) 
iterative solution of the DMFT equations on the imaginary axis using QMC, iii) computa- 
tion of observables of interest. This procedure has been repeated for several lattice parame- 
ters corresponding to pressures up to several 100 GPa. In the case of MnO we employed the 
Hirsch-Fye QMC algorithm available at the time, for which the calculations were limited to 
the temperature T = 1160 K. For the Fe203 study the continuous-time QMC algorithm (hy- 
bridization expansion) was used, which allowed us to reach room temperature and study the 
anti-ferromagnetic ordering. In both cases we had to introduce special Monte-Carlo moves to 
ensure ergodic simulation for parameters close to the spin transition. The calculated observ- 
ables include the single-particle spectral density to study the formation of the charge gap, 
local orbital occupancies (rij) and static local density-density correlations {niUj) to investigate 
the orbital redistributions, instantaneous local moments and interaction energy, and local spin 
susceptibility to determine the screening of the local spin moments. 

In Fig.|4]we show the evolution of the local moments and the orbital occupancies as functions 
of specific volume. We use two different quantities to characterize the local moment: (a) the 
mean instantaneous moment defined as an equal time correlation function Ms = \J (™^) and 
(b) the screened local moment defined through the local spin susceptibility Moff = -y/Txioc- 
These two different definitions of the local moment yield almost the same T-independent values 
in materials with Curie- Weiss behavior. The two studied materials exhibit clear similarities. In 
particular, the HS state with Mg ^ 5 and orbital occupancies corresponding to the complete spin 
polarization of a half-filled d shell dictated by the first Hund's rule, the LS state characterized 
by emptying of the eg orbitals, substantial reduction of the local moment and the screening of 
the local moment refiected in the difference between Mg and A/off. Furthermore it turns out 
that the charge gap disappears in the LS phase (see Figs. [5] and [51). 

Nevertheless, there are also several differences between MnO and Fe2 03. Their orbital oc- 
cupancies differ quantitatively, especially the Cg occupancies in the LS phase are substantially 
larger in Fe203. This is a result of more covalent character of this compound, i.e., stronger 
p-Cg mixing. Another consequence of the strong covalent p-Cg bonding is a transfer of charge 
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Fig. 5. The single-particle spectra of Fe203 at various specific volumes (T — 580 K) [44]. The HS 
solutions are shown in the left and LS solutions in the right panel. The Fe d spectra are resolved into 
the eg (called eg in the legend) and t^g (further split in aig and eg by distortion from precise octahedral 
symmetry) contributions. The O p spectra are marked by blue shading. 




between Fe and O upon the HS-LS transition, reflected in the difference of the total occupancies 
Ud in the two phases. There is also a qualitative difference between the transitions in MnO and 
Fe2 03. In MnO we find a smooth crossover characterized by a unique solution for all specific 
volumes with intermediate magnetic moments, and only when the lattice is included we observe 
a phase separation into high- volume HS phase and low- volume LS phase in the crossover region. 
In Fe203 we observe an electronically first-order transition with hysteresis. 

Next, we compare the evolution of the single-particle spectral densities (see Figs. [S] and IH]) as 
the volume is reduced. In both cases we observe the closing of the HS charge gap as the material 
is compressed, and LS spectra that bear clear resemblance to the non-interacting bands with 
some quasiparticle renormalization. However, there is a clear difference in the way the HS gap 
disappears. In Fe203 (Fig. [5]) the gap can be squeezed arbitrarily small without qualitative 
changes of the spectra and the HS phase becomes unstable only when the gap is zero (within 
the numerical accuracy and resolution). In contrast, the transition in MnO does not proceed by 
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Fig. 7. Comparison of the phase boundaries of the two-band obtained with full DMFT (QMC) calcu- 
lation (taken from Ref. [55]) with the phase boundaries obtained by simple ionic estimates described in 
the text. The boundary between the "metal" and "band insulator" was obtained by setting the charge 
gap above the band insulator ground state to zero. 

squeezing the gap to zero, but by appearance of in-gap states. We use this observation below 
to discuss how the HS phase is destabilized. 

For the sake of clarity we simplify the local interaction to ^ niuinj-^ ~\- [U — J)nianja- 
in the following discussion. We will estimate two quantities, the charge gap in the HS phase: 

Eg « E{d^) + E{d^) ~ 2E{d'^) ~W = U + 4J - A-W, (3) 

and the difference (per atom) between the LS and HS ground states: 

£;(LS) - £;(HS) = 6 J - 2Z\. (4) 

Here A is the eg-t2g crystal-field splitting, W is an effective bandwidth, which describes reduc- 
tion of the ionic gap due to the actual bandwidth and the charge-transfer parameter, E{d^) 
is the ground-state energy of an ion in the HS state and E{d^), E{d'^) are the energies of the 
corresponding electron addition and removal states, respectively. The estimates are obtained 
by simple ionic considerations. The interaction parameters U and J are considered to remain 
constant, while the crystal-field splitting A and the bandwidth W increase with pressure. The 
insulating HS state is destabilized when either ([3]) or (|4]) becomes zero. Which of the two quan- 
tities reaches zero first depends on the parameters of the system. We have adopted the terms 
gap closing and local state transition for the former and the latter instability, respectively. 

Comparing the behavior of the charge gap in Fe203, where the HS phase exists as long as 
the gap is finite, and in MnO, where the transition starts while the gap is still finite, leads 
to the conclusion that the observed metal-insulator transition in Fe203 proceeds by the gap 
closing while in MnO the local state transition takes place. Interestingly, a similar behavior 
was observed in the two-band model studied by Werner and Millis . In Fig. [7] we put the 
transition lines obtained from approximate expressions like ^ and ^ adapted for the two- 
band model on top of the numerical results obtained by the DMFT solution of the two-band 
model. Thus we can identify the transition between the "Mott insulator" and the "metal" , in 
the language of Ref. [92], with the gap closing and the transition between "Mott insulator" 
and the "band insulator" with the local state transition. Note that the third transition (black 
line) exists in the two-band model with two electrons, but not in a five-band model with five 
electrons relevant for the studied materials. The example of the two-band model shows that the 
gap closing and local state transition concepts are quite general. Similar ionic considerations 
were recently used by Lyubutin et al. [74 to other fillings of the d shell. It was also shown that 
for d^ Eq. [3| holds only if Z\ < 2 J. 

In studying the relationship between metal-insulator and spin transitions in Fe203 and 
MnO we discussed two scenarios how an insulating HS phase can be destabilized. In the first 
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one, the charge gap is closed due to increase of the crystal field as well as the bandwidths so 
that the system becomes metallic, at which point the ionic reasoning presented here and in 
Ref. 74; breaks down. In the second one, the local HS becomes energetically unfavorable while 
the charge gap is still finite and the system transforms into a LS phase, which may be either 
metallic or insulating, e.g., LuFeOs |93l94j . depending on the parameters of the system. 

4 Correlation-induced atomic displacements and structural transformations 

Previous applications of the LDA+DMFT approach have focused on the investigation of cor- 
relation effects in the electronic system with a given lattice structure, thereby neglecting the 
mutual interaction between electrons and ions. As a result, the influence of the electrons on the 
lattice structure was disregarded. Recent applications of LDA+DMFT, e.g., computations of 
the volume collapse in paramagnetic Ce [30 31 32, and Pu |33I34I35] and of the collapse of the 
magnetic moment in MnO 45', did include the lattice, but only calculated the total energy of 
the correlated material as a function of the unit cell volume [53]. Moreover, these implementa- 
tions of LDA-I-DMFT cannot describe the electronic and structural properties of electronically 
correlated materials in the case of structural transformations, e.g., involving the cooperative 
Jahn- Teller (JT) effect [96 97 98 99 100 . This is due to the atomic-sphere approximation within 
the linearized and higher-order muffin-tin orbital [L(N)MTO] techniques [1011102) which are 
conventionally used in the LDA-I-DMFT approach. The approximation employs a spherical po- 
tential inside the atomic sphere, thereby neglecting multipole contributions to the electrostatic 
energy due to the distorted charge-density distribution. This makes the L(N)MTO technique 
unsuitable to determine atomic positions reliably. Instead, the recently proposed implementa- 
tion of the LDA+DMFT approach [103I104I105I106I107I108] . which employs plane- wave pseu- 
dopotentials [10111] and thus avoids the atomic sphere approximation, does not neglect such 
contributions. Thereby it becomes possible to describe the effect of the distortion on the elec- 
trostatic energy [1031104] . 

Here we discuss the LDA+DMFT scheme implemented with plane-wave pseudopoten- 
tials [10311041105] to compute structural transformations, e.g., structural phase stability and 
structure optimization, caused by electronic correlations. Most importantly, this scheme is able 
to determine correlation-induced structural transformations in both paramagnetic and long- 
range ordered solids. Thereby it is able to overcome the limitations of standard band-structure 
approaches and earlier implementations of the LDA+DMFT approach, and opens the way for 
fully microscopic investigations of the structural properties of strongly correlated electron mate- 
rials. In the following we present the application of this approach to the investigation of orbital 
order and the cooperative JT distortion in two prototypical JT materials, KCuFs and LaMnOa, 
and compute the electronic, structural, and orbital properties in their room-temperature para- 
magnetic phase. These are the first results obtained for a structural optimization where the 
stability of the cooperative JT distortion in paramagnetic KCuFa and LaMnOa is investigated 
using total energy calculations jl09] . 

4.1 Application to KCuFs 

KCuFa is a prototype material regarding the cooperative Jahn- Teller (JT) effect, orbital or- 
dering, and low-dimensional magnetism [97I98I99I100I110] . It is an insulating pseudocubic per- 
ovskite whose structure is related to that of high-temperature superconductors and colossal 
magnetoresistance manganites. The copper ions have octahedral fluorine surrounding and are 
nominally in a Cu^+ (3d^) electronic configuration, with completely filled t2g orbitals and a sin- 
gle hole in the Cg states. This gives rise to a strong JT instability that lifts the cubic degeneracy 
of the Cu eg states due to a cooperative JT distortion [97I98I99I100] . The latter is characterized 
by CuFg octahedra elongated along the a and b axis and arranged in an antiferro-distortive 
pattern in the ab plane [111] . The strong JT distortion persists up to the melting temperature 
(~ 1000 K) and is associated with the particular orbital order in KCuFa, in which a single hole 
alternatingly occupies d^^-z^ and dy2_z2 orbital states, resulting in a tetragonal compression 
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(c < a) of the unit cell. Purely electronic effects as in the Kugel-Khomskii theory |97l98l99ll00j 
and the electron-lattice [112' interaction have been discussed as a possible mechanism behind 
the orbital ordering in KCUF3. Nevertheless, the mechanism responsible for the orbital order 
in KCuFg is still being debated in the literature [8l97l98l99ll()()ll():^lll2lll8lll4lll,^j . 

KCuFs has a relatively high (tetragonal, space group I A/ mem) crystal symmetry, hence it 
is one of the simplest systems to study. The JT distortion in KCuFa can be expressed using 
a single internal structure parameter, the shift of the in-plane fluorine atom from the Cu- 
Cu bond center. Moreover, KCuFa has a single hole in the 3d shell resulting in the absence 
of multiplet effects. The electronic and structural properties of KCuFs have been intensively 
studied by density functional theory in the local density approximation (LDA) [3] , the general- 
ized gradient approximation (GGA) [5161116] . or using the so-called LDA-I-U approach |7I8I9| . 
While the LDA-f U calculations account rather well for the value of the equilibrium JT dis- 
tortion in KCuFa |115) , the calculations simultaneously predict a long-range antiferromagnetic 
order which indeed occurs in KCuFa below Tn ~ 22-38 K [117j . The LDA-I-U calculations 
give the correct insulating ground state with the long-range A-type antiferromagnetic and 
dx^-z^ / dy2-z^ antiferro-orbital order |8lll3llf5] . consistent with the Goodenough-Kanamori- 
Anderson rules for a superexchange interaction. Moreover, LDA-I-U calculations for a model 
structure of KCuFs in which cooperative JT distortions are completely neglected reproduce the 
correct orbital order, suggesting an electronic origin of the ordering |8lll3j in agreement with 
the Kugel-Khomskii theory [97 98 99 100 . Nonmagnetic LDA/GGA calculations instead pre- 
dict a metallic behavior and cannot explain the insulating paramagnetic behavior at T > Tn. 
The electronic and structural properties of KCuFa have been recently reexamined by means of 
LDA-I-U molecular-dynamic simulations, indicating a possible symmetry change and challenging 
the original assignment of tetragonal symmetry |115| . This symmetry change seems to allow 
for a better understanding of Raman [118j . electronic paramagnetic resonance [11911201121) . 
and x-ray resonant scattering [1221123] properties at T w Tat. However, the details of this dis- 
tortion have not been fully understood yet. While the LDA-I-U approach is able to determine 
electronic properties and the JT distortion in KCuFs rather well [115j . its application is lim- 
ited to temperatures below Tn- Therefore LDA-I-U cannot explain the properties at T > Tn 
and, in particular, at room temperature, where KCuFa is a correlated paramagnetic insulator 
with a robust JT distortion which persists up to the melting temperature. To determine the 
correct orbital order and cooperative JT distortion for a correlated paramagnet, i.e., to perform 
a structural optimization, we here employ the novel GGA-I-DMFT scheme implemented with 
plane- wave pseudopotentials [10311041105] . 

We start by calculating the nonmagnetic GGA band structure of KCuFa within the plane- 
wave pseudopotential approach [124] . In this calculation we use the Perdew-Burke-Ernzerhof 
exchange-correlation functional together with Vanderbilt ultrasoft pseudopotentials for copper 
and fluorine, and a soft Troullier-Martin pseudopotential for potassium. The calculation was 
performed for the experimental room-temperature crystal structure with space group I A/ mem 
and lattice constants a = 5.855 and c — 7.852 A [111] . We have used different values of the JT 
distortion 5jt defined accordingly as (5jt = \{di — dg) / {di + ds)\ di and dg denote the long and 
short Cu-F bond distances in the ah plane of CuFe octahedra, respectively, and 2{di + dg) — a. 
We express the distortion Jjt in percent of the lattice constant a, e.g., Sjt — 0.002 = 0.2%, and 
consider 0.2% < Jjt < 7%. The structural data [111] at room-temperature yield (5jt = 4.4%. In 
the present calculation we keep the lattice parameters a and c and the space group symmetry 
fixed. 

For all values of 5jt considered here, the nonmagnetic GGA calculations yield a metallic 
rather than the experimentally observed insulating behavior. We also find an appreciable orbital 
polarization due to the crystal field splitting. Overall, the GGA results qualitatively agree with 
previous band-structure calculations [81115] . Thus, the electron-lattice interaction alone is found 
insufficient to stabilize the orbitally ordered insulating state. The GGA total energy profile is 
seen to be almost constant for 6jt < 4%, with a very shallow minimum at about 2.6%. (see 
Fig. [51 left panel). This would imply that KCuFa has no JT distortion for temperatures above 
100 K, which is in clear contradiction to experiment. Obviously, a JT distortion by itself. 
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Fig. 8. Left panel: Comparison of the total energies of paramagnetic KCuFs computed by GGA and 
GGA+DMFT (QMC) as a function of the JT distortion [1031104) . Error bars indicate the statistical 
error of the DMFT(QMC) calculations. Right panel: Orbitally resolved Cu eg spectral densities of 
paramagnetic KCuFs as obtained by GGA+DMFT(QMC) for different values of the JT distortion. The 
resulting orbitally resolved spectral density which is shown here by solid [dashed] line is predominantly 
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character (in the local frame |125p . 



without the inclusion of electronic correlations in the paramagnetic phase, cannot explain the 
experimentally observed orbitally ordered insulating state of KCuFs. 

The next step is the construction of an effective lovsr-energy Hamiltonian Hqqa for the 
correlated, partially filled Cu eg orbitals for each value of the distortion i5jt considered here. This 
is achieved by employing the pseudopotential plane- wave GGA results and making a projection 
onto atomic-centered symmetry-constrained Cu eg Wannier orbitals [511521103] . Taking the 
local Coulomb repulsion U = 7 eV and Hund's rule exchange J = 0.9 eV into account, we 
obtain the following low-energy Hamiltonian for the two (m = 1,2) Cu eg bands: 

H = Hgga + U'^ fiimtnimi + ^ 5„„,J)nii„hi2c,' - Hue- (5) 

1771 iaa' 

Here the second and third terms on the right-hand side describe the local Coulomb interaction 
between Cu eg electrons in the same and in different orbitals, respectively, with V — U — 2 J, 
and i?DC is a double counting correction which accounts for the electronic interactions al- 
ready described by the GGA (see below). We solve the many-body Hamiltonian ^ for each 
value of JjT using the single-site DMFT with Hirsch-Fye quantum Monte Carlo (QMC) cal- 
culations (53il26| . The calculations were performed at T = 1160 K (/3 = 10 eV~^), using 40 
imaginary-time slices |127j . For all values of Sjt considered here, we compute the GGA+DMFT 
total energy as [3211031104] 

E = Egga [p] + (Hgga) - J] e^,^,^ + (Hu) - Euc , (6) 

Tn.k 

where Eqqa[p] is the total energy obtained by GGA. The third term on the right-hand side 
of Eq. (jnj) is the sum of the GGA Cu eg valence-state eigenvalues and is given by the thermal 
average of the GGA Hamiltonian with the GGA Green function G^'^^(ia;„): 

E ''^T = ^ E Tr[i^GGA(fc)G^^(^.„)]e'-°^ (7) 

771. k 71. k 

The average {Hqqa} is evaluated similarly but with the full Green function including the 
self-energy. The interaction energy (Hu) is computed from the double occupancy matrix. The 
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Fig. 9. The M/mcm primitive cell and hole orbital order as obtained by the GGA+DMFT calcula- 
tion for paramagnetic KCuFa with Sjt = 4.4% is shown. The fluorine atoms and fluorine octahedra 
are shown in green, the potassium in red, and the Wannier Cu eg charge density in blue. The local 
coordinate system is chosen with the z direction deflned along the longest Cu-F bond of the CuFe 
octahedron. 

double-counting correction Euc = \U N^^iN^^ — 1) — ^JN^giNe^ — 2) corresponds to the 
average Coulomb repulsion between the N^^ electrons in the Cu Cg Wannier orbitals. Since the 
Hamiltonian involves only correlated orbitals the number of Wannier electrons Nd is conserved. 
Therefore, the double-counting correction amounts to an irrelevant shift of the total energy. 

In Fig. [5] (left panel) we present our results for the GGA-f-DMFT total energy as a function 
of the JT distortion Jjt- The inclusion of the electronic correlations among the partially filled 
Cu eg states in the GGA-f DMFT approach is seen to lead to a very substantial lowering of the 
total energy by ~ 175 meV per formula unit. This implies that the strong JT distortion persists 
up to the melting temperature (> 1000 K), in agreement with experiment. The minimum of 
the GGA-I-DMFT total energy is located at the value Sjt = 4.2%, which is also in excellent 
agreement with the experimental value of 4.4% [111) . This clearly shows that the JT distortion 
in paramagnetic KCuFs is caused by electronic correlations. 

In Fig. \E\ (right panel) we show the spectral density of paramagnetic KCuFa calculated for 
several 5jt values using the maximum entropy method. The GGA-f-DMFT calculations give 
a paramagnetic insulating state with a substantial orbital polarization for all (5jt considered 
here. The energy gap is in the range 1.5-3.5 eV, and increases with increasing of (5jt- The 
sharp feature in the spectral density at about —3 eV corresponds to the fully occupied 3z^ — 

orbital [125] . whereas the lower and upper Hubbard bands are predominantly of — y'^ 
character and are located at —5.5 eV and 1.8 eV, respectively. The corresponding Cu Cg Wannier 
charge density calculated for the experimental value of JT distortion of 4.4% is presented in 
Fig. [21 The GGA-f DMFT results clearly show an alternating occupation of the Cu d^^-za 
and dy2_z2 hole orbitals which implies antiferro-orbital order. Moreover, we have proved the 
stability of the paramagnetic solution at high temperatures (560 K) with respect to the A-type 
antiferromagnetic one. In this calculation we have used the experimental room-temperature 
crystal structure of KCuFa with (5jt — 4:A% and long-range A-type antiferromagnetic order. 
However, in agreement with experiment, the calculation gives paramagnetic insulating solution 
with the orbital order as it has been found above. 

We now perform a structural optimization of paramagnetic KCuFa. For simplicity, we con- 
sider only two independent structural parameters, the lattice constant a and the JT distortion 
(5jT. In this calculation we keep unchanged the space group symmetry (tetragonal 1 4/ mem) 
and the experimental value of the unit cell volume (taken at the ambient pressure at room 
temperature) [1 IJJ . We first calculate the non- magnetic GGA electronic structure for different 
values of 6jt and lattice constant a. In Fig. [TU] (left panel) we display the GGA total energies 
for different JT distortion 6jt as a function of the lattice constant a. The data points were 
further interpolated by smooth curves. The result of the total energy variation, namely, the line 
connecting the minima of the corresponding curves, is shown in bold (blue) line. We find that 
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Fig. 10. Comparison of the total energies of paramagnetic KCuFa computed by GGA (left panel) and 
GGA+DMFT(QMC) (right panel) for different values of the JT distortion Sjt as a function of the 
lattice constant a [1U4) . The result of the total energy variation is shown as the thick blue line. Error 
bars indicate the statistical error of the DMFT(QMC) calculations. 

the variation of the GGA total energy, is seen to be constant for (5jt < 2% with the end point 
at a ~ 5.75 A. This imphes the absence of the cooperative JT distortion and results in a nearly 
cubic (c/a ~ 1.0) unit cell, which is in clear contradiction to experiment [lllj . 

To proceed further, we construct the effective low-energy Hamiltonian for the partially filled 
Cu eg orbitals for each value of the JT distortion Sjt and the lattice constant a considered here, 
and compute the corresponding GGA+DMFT total energy (see Fig. [TUJ right panel) |103I104| . 
In contrast to the GGA result, the inclusion of the electronic correlations among the partially 
filled Cu Cg states in the GGA+DMFT method not only correctly describes the spectral prop- 
erties, but also leads to a very prominent minimum in the resulting total energy variation. 
The minimum is located at the value a = 5.842 A and Sjt « 4.13%, which is in excellent 
agreement with experimental value a = 5.855 A and Sjt ~ 4.4%. In addition, the structural 
optimization within GGA+DMFT predicts the correct tetragonal compression of the unit cell 
c/a « 0.95 [TTT] . 



4.2 Application to LaMnOs 

LaMnOs is another prototypical material with respect to the cooperative Jahn- Teller effect and 
orbital order. At room-temperature it is an orthorhombic perovskite with a GdFeOs-like crystal 
structure and space group Pnma |128] . The Mn ions have octahedral oxygen surrounding and 
are in a high-spin 3d'' electronic configuration due to Hund's rule coupling, with three electrons 
in the t2gf orbitals and a single electron in an Cg^ orbital (^ig^g orbital configuration). The 
two types of structural instabilities give rise to the changes relative to the cubic perovskite 
structure. The first is related to a large ion-size misfit parameter -j2{Ro + RMn)/{Ro + Rlo) 
which favors rotations of the octahedra to accommodate a more efficient unit cell space filling; 
Rivim Rha, and Ro are the ionic radii of Mn, La, and O ions, respectively. The second is a JT 
instability due to the e^-orbital degeneracy which is lifted due to a cooperative JT distortion, 
leading to an alternating occupation of d^^2_j.2 and d^y2_^2 electron orbitals in the ab plane 
(antiferro-orbital ordering) and to a tetragonal compression of the unit cell. The rotations of the 
octahedra lower the symmetry further, finally resulting in the orthorhombic unit cell. The JT 
distortion experimentally persists up to Tjt ~ 750 K which is remarkably higher than the Neel 
temperature of ~ 140 K. Above Tjt LaMnOa undergoes a structural phase transition jl29) . 
with volume collapse |130) to a nearly cubic structure in which orbital order and JT distortion 
vanish [129) . Below T/v ~ 140 K LaMnOa shows A-type long-range antiferromagnetic order 
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consistent with the Goodenough-Kanamori- Anderson rules for a superexchange interaction with 
dy,,j;2_^2 / dj,y2_j.2 antiferro-orbital order |128ll29j . 

The LDA+U approach is found to give, at equihbrium, the correct insulating behavior 
of the low-temperature antiferromagnetic phase and a JT distortion in satisfactory agree- 
ment with experiment (T?l j . Nevertheless, it should be noted again that the application of 
this approach is limited to temperatures below Tjy. Therefore, LDA-I-U cannot describe the 
properties of LaMnOs at T > Tjsi and, in particular, at room temperature, where LaMnOs 
is a correlated paramagnetic insulator with a robust JT distortion. Nonmagnetic LDA/GGA 
calculations instead give a metal and cannot explain the insulating paramagnetic behavior at 
T > Tjv- The electronic properties of paramagnetic LaMnOa have already been studied using the 
LDA-I-DMFT approach 132 133 134 135 . However, no attempt has been made to determine the 
structural properties and, in particular, the value of the cooperative JT distortion so far. Here 
we investigate the stability of the cooperative JT distortion and perform a structural optimiza- 
tion of paramagnetic LaMnOa. To determine the electronic and structural properties we employ 
the novel GGA+DMFT scheme implemented with plane- wave pseudopotentials |103I104I105) . 

We first compute the nonmagnetic GGA band structure of paramagnetic LaMnOs using the 
plane- wave pseudopotential approach as in the case of KCuFs. In this calculation, we have used 
the orthorhombic Pnma crystal structure as reported by Elemans et al, with lattice constants 
a = 5.742, b = 7.668, and c = 5.532 A |128j . The calculation was performed for different values 
of JT distortion (5jt, which is now defined as the ratio between the difference of the long (di) 
and the short (ds) bond distances and the mean Mn-0 distance in the basal ab plane, i.e., 
(5jT — ^{di — ds)/{di + ds). Structural data |128) yield (5jt = 0.138. Here we change only the 
parameter Jjt (0 < (5jt < 0.2) while the value of the MnOg octahedron tilting and rotation 
was fixed. 

For all values of (5jt considered here we find a metallic solution with a considerable orbital 
polarization due to the crystal field splitting. Overall, these results qualitatively agree with 
previous band structure calculations |136) . namely, GGA cannot describe a paramagnetic insu- 
lating behavior which is found in experiment. We notice that even for the large 5jt value of 0.2 

45% larger than found in experiment [128] ) the GGA calculations predict a metal. In Fig.fTTI 
(left panel) we display our results for the GGA total energy as a function of the JT distortion 
(5jT. The GGA total energy is almost parabolic which implies the absence of a cooperative 
JT distortion. Altogether, the non-magnetic GGA calculations give a metallic solution without 
cooperative JT distortion that contradicts experiment [1281129] . This implies the importance of 
electronic correlations to describe the experimentally observed insulating behavior and orbital 
order in paramagnetic LaMnOa. 

To proceed further, we turn to the GGA-I-DMFT results where we treat the Mn eg orbitals 
as correlated orbitals. To achieve this we employ the GGA results and make a projection onto 
atomic-centered symmetry-constrained Mn Cg Wannier orbitals |51I52I105] . In this calculation 
we assume that three (among the 3d^ electronic configuration) electrons are localized in the 
t2g orbitals. Therefore, they are treated as classical spins S*, with a random orientation above 
Tjv (i.e., there is no correlation between different Mn sites), which couple to the Cg electron 
with an energy JS. This coupling can be estimated as the energy of the splitting of the eg-|- and 
Bgj, bands in the ferromagnetic band-structure calculations and gives an additional term in the 
Hamiltonian ([S]), namely, 

H = HqGA + "-""f^jniJ- ~ JS'^^flira^ ~ fliral) + '^(y - 6aa' J)nilant2a' ~ HuC ■ (8) 

im im iaa' 

This corresponds to the ferromagnetic Kondo-lattice model Hamiltonian with an on-site 
Coulomb repulsion between eg electrons, which has been intensively studied as a possible mi- 
croscopic model to explain the colossal magnetoresistance in manganites |138I139) . We note 
that in order to calculate the total energy one needs to modify Eq. [Blby adding the expectation 
value of the JS term which describes the total energy gain due to the spin polarization of the 
eg orbitals at the Mn site. We take the local Coulomb repulsion [/ = 5 eV, the Hund's rule 
exchange J = 0.75 eV, and 2JS = 2.7 eV from the literature (1331134] . and further solve the 
many-body Hamiltonian for each value of Sjt using the single-site DMFT with Hirsch-Fye 
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Fig. 11. Left panel: Comparison of the total energies of paramagnetic LaMnOs computed by GGA 
and GGA+DMFT (QMC) as a function of the JT distortion |104j . Error bars indicate the statistical 
error of the DMFT(QMC) calculations. Right panel: Orbitally resolved Mn eg spectral densities of 
paramagnetic LaMnOs as obtained by GGA+DMFT(QMC) for different values of the JT distortion. 
The resulting orbitally resolved spectral density shown by solid [dashed] line is predominantly of 3z^ — r^ 
[x^ — y^] character (in the local frame [IS?) ). 



quantum Monte Carlo (QMC) calculations [531126) . Again the calculations were performed at 
T = 1160 K (/3 = 10 eV"^), using 40 imaginary-time slices |127) . 

In Fig. [11] (left panel) we display the GGA-I-DMFT total energy calculated as a function 
of the JT distortion Sjt- In contrast to the GGA result, GGA+DMFT yields a substantial 
total energy gain of ~ 150 meV per formula unit. This implies that the cooperative JT dis- 
tortion indeed persists up to high temperatures (T > 1000 K), while in GGA a JT distortion 
does not occur at all. Taking into account that the calculations have been performed for the 
low-temperature crystal structure of LaMnOa [128j this estimate (150 meV) is in good agree- 
ment with TjT ^ 750 K at which the JT distortion vanishes [129I130[ . The minimum of the 
GGA+DMFT total energy is located at the value Sjt ~ 0.11, which is also in good agreement 
with the experimental value of 0.138 [128I129| . We note that GGA+DMFT calculations cor- 
rectly describe both electronic and structural properties of paramagnetic LaMnOa. This shows 
that the JT distortion in paramagnetic LaMnOa is caused by electronic correlations. 

In Fig. [TT] (right panel) we display the spectral density of paramagnetic LaMnOa calculated 
for several S]t values using the maximum entropy method. We find a strong orbital polarization 
for large (5jt, which gradually decreases for decreasing JT distortion. The occupied part of the 
eg density is located at about —2 to —1 cV and corresponds to the Cg states with spin parallel 
to the t2g spin at that site. It has predominantly Mn d^^2_j,2 and d^y2_j.2 orbital character with 
a considerable admixture of d^2_^2 for small JT distortions. The energy gap is about 2 eV for 
large Sjt and considerably decreases with decreasing Sjt, resulting in a pseudogap behavior at 
the Fermi level for Sjt = 0. In Fig. [T^ we present the corresponding Mn eg Wannier charge 
density computed for the experimental JT distortion value of Sjt = 0.138. The result clearly 
indicates an alternating occupation of the Mn (1^x2 ^^2 and d^y2_j.2 orbitals, corresponding to 
the occupation of a iz^ — orbital in the local frame [137[ . which implies antiferro-orbital 
order. Thus, in agreement with experiment, the calculations give a paramagnetic insulating 
solution with antiferro-orbital order and stable JT distortion. 

To conclude, we presented applications of the novel GGA+DMFT approach to two pro- 
totypical Jahn- Teller materials, KCuFa and LaMnOa, and computed the orbital order and 
cooperative JT distortion in these compounds. In particular, our results obtained for the para- 
magnetic phase of KCuFa and LaMnOa, namely an equilibrium Jahn- Teller distortion Sjt of 
4.2% and 0.11, respectively, and antiferro-orbital order, agree well with experiment. The present 
approach overcomes the limitations of the LDA+U method and previous implementations of 
LDA+DMFT and is able to determine correlation-induced structural transformations in both 
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Fig. 12. The Pnma primitive cell and orbital order as obtained by the GGA+DMFT calculation for 
paramagnetic LaMnOs with Sjt ~ 0.138 is shown. The oxygen atoms and oxygen octahedra are shown 
in green, the lanthanum in red, and the Wannier Mn eg charge density in blue. 

paramagnetic and long-range magnetically ordered solids. Thereby it can be employed for the 
lattice optimization and molecular-dynamic simulations of these systems. The computational 
scheme presented here allows us to investigate complex materials with strong electronic interac- 
tions and correlation-induced structural transformations, shifts of equilibrium atomic positions 
and changes in the lattice structure, and to perform a structural optimization of paramag- 
netic solids. It opens the way for fully microscopic investigations of the structural properties 
of strongly correlated electron materials such as lattice instabilities observed at correlation- 
induced metal-insulator transitions. 



5 Kinks in the dispersion of correlated materials 

The electronic dispersion relation indicates at which energy and crystal momentum one-particle 
excitations can occur in a solid. Interactions between electrons influence not only the shape and 
position of the dispersion, but also make the lifetime of such excitations finite. A robust general 
definition for the branches of the dispersion relation is provided by the local maxima of 
the spectral function A{k,u!), i.e., its peak locations as a function of uj for a given momentum 
k. For energies w far from the Fermi surface, the peaks in the spectral function of a strongly 
interacting material can be quite broad, corresponding to incoherent excitations with a short 
lifetime. Nevertheless the maxima in A{k,ijj) are always well defined, even if their location might 
be difficult to determine from angle-resolved photoemission data at high binding energies. The 
excitations for energies close to the Fermi surface, on the other hand, are coherent in a Landau- 
Fermi liquid, i.e., they are only weakly damped and correspond to sharp peaks in the spectral 
function. In this case the correlated dispersion is given by the uncorrelated dispersion multiplied 
by a renormalization factor Z related to Fermi-liquid parameters. 

A characteristic feature of strongly correlated metals is a transfer of spectral weight into 
the lower and upper Hubbard bands due to the Coulomb interaction. Together with the central 
peak near the Fermi energy the momentum-integrated spectral function therefore typically has 
three peaks. One exemplary material with this type of spectrum is the cubic perovskite SrVOs 
with 3c?^ electronic configuration, leading to band-filling of one electron per site in the threefold- 
degenerate t2g bands. To obtain its momentum-resolved spectral function with LDA-f DMFT, 
we used iVth order muffin orbitals (NMTO) to obtain the downfolded 3x3 Hamiltonian for the 
t2g bands [WO,. For the DMFT calculation the Coulomb interaction parameters for these bands 
are U = 5.55 eV for the on-site repulsion, J = 1.0 eV for the Hund's exchange coupling, and U' = 
3.55 eV for the interorbital Coulomb repulsion |141) . The local (momentum-integrated) spectral 
function has an asymmetric three-peak structure with a high central peak and two smaller peaks 
due to the Hubbard bands [140j . The total spectral function A{k,uj) is shown in the intensity 
plot in Fig. 1131 with its local maxima marked by white dots. Remarkably, the dispersion Ek 
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Fig. 13. Intensity plot of the total spectral function for SrVOs obtained with LDA+DMFT; after [HD]. 
Near the Fermi energy the dispersion Ek (white dots) is given by the renormalized LDA band structure 
Efe (black lines) with a Fermi-liquid factor Z = 0.35, i.e. Ek = Zek (blue line). For higher energies the 
dispersion relation follows the LDA band structure with a different renormalization, Ek = Z'ek + c± 
(pink line), with Z' = 0.64, c+ = 0.086 eV, c_ = 0.13 eV |142j. The crossover between the two slopes 
leads to kinks at energies 0;*,+ = 0.22 eV and tu*,- = -0.24 eV in the effective dispersion. The parameter 
values given here are calculated from the theory of Ref. |142) : due to the shift c± the parameter Z' 
differs from the fitted parameter in Ref. |140j without shift. 



shows a small sudden change in slope near energies uji, ~ 0.25 eV. There a crossover between 
two regimes occurs: In the low-energy Fermi-liquid regime, \uj\ < cj^,, the LDA dispersion ek 
is renormalized by the Landau parameter Z so that Ef. = Ztk- On the other hand, at higher 
energies, outside the Fermi-liquid regime but still inside the central peak, the dispersion is 
given by = Z'e^ with a different renormalization factor. In Ref. [140) these parameters were 
determined as Z — 0.33 and Z' = 0.53, i.e., the renormalization is much stronger near the Fermi 
surface. This value of Z' w 0.5 was also found in ARPES experiments [ 143.144 . Although no 
sharp kink was found, a weak change in slope was detected at w w —0.2 eV |144j . in agreement 
with the LDA-FDMFT calculation. 

These kinks in the dispersion can be understood from a microscopic theory |142] . which 
correctly predicts the renormalized bandstructure outside the Fermi-liquid regime only from 
Fermi-liquid parameters and the uncorrelated spectrum. A sufficiently strong interaction is 
needed, so that the local spectral function has a three-peak structure with a central peak and 
two Hubbard peaks. In the absence of particle-hole symmetry this condition can be relaxed; 
e.g., it suffices that the central peak is well separated from the lower Hubbard peak for a 
kink to develop below the Fermi surface. As we will now discuss, the dips on the sides of the 
central peak lead to kinks in the dispersion inside the central peak; for simplicity, we consider 
only a one-band system with symmetric spectrum (Fig. [H)). Namely, well-separated peaks in 
the local spectral function A{u!) = — Im[G(LLi)/7r] (top panel in Fig. [H)) have the consequence 
that due to Kramers-Kronig relations, Re[G'(w)] must have sharp maxima at energies Wmax 
inside the central spectral peak (second panel from top in Fig. ll4p . These maxima induce kinks 
in the self-energy E{k,uj), which lead to kinks in the fe-resolved spectral function, A{k,uj) 
= — Im[w + n — ek — S{k,uj)]~^ /tt, and thus in the dispersion. In DMFT the self-energy is 
independent of momentum, S{k, lu) — and related to the local Green function by the self- 

consistency condition, X!{uj) = to + ^ — 1/G{uj) — A{G{u!)), where A{G{uj)) is the hybridization 
function. For a three-peak spectral function A{uj), the term Re[w + /i — 1/G{ijj)] is essentially 
linear inside the central peak (bottom panel in Fig. [14]). By contrast, A{G{u!)), which can be 
expanded in powers of G{ui), always contains a term linear in G{uj). As a consequence S{io) 
always has a contribution from the peaks in Re[G'(a;)] that is nonlinear in oj inside the central 
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Fig. 14. Origin of kinks in a strongly correlated system (one-band Hubbard model in DMFT with 
Z — 0.086) [142| . From top to bottom: Local spectral function with three well-separated peaks; corre- 
sponding real part of the local Green function with peaks inside central spectral peak; real part of the 
self-energy, in which the maxima of Re[G(cj)] translate into kinks; uj — 1/G(a;), which contributes only 
a smooth linear background to Re[X'(aj)]. 



peak. Thus Re[Z'(a;)] has two different slopes for \Ld\ smaller or larger than UJ^,. This translates 
into the usual Fermi-liquid renornialization = Zef. near the Fermi energy, but a shifted 
and renormalized dispersion = Z'e^ + const for higher binding energies |142j . From the 
Fermi-liquid renormalization Z and the non-interacting bandstructure the constant shift, the 
parameter Z' and the crossover energy can be obtained explicitly |142| ; in particular Z' and 
are proportional to the Fermi-liquid parameter Z (multiplied by an uncorrelated scale) . The 
values for SrVOs are listed in Fig. [T3] and fit the correlated dispersion very well. 

According to the theory described above the kinks in the slope of the dispersion are a direct 
consequence of the electronic interaction [142) . They can lead to corresponding kinks in the 
low-temperature electronic specific heat [145) . The presence of kinks has also been linked to 
maxima in the spin susceptibility |146| . Of course, additional kinks in the electronic dispersion 
may arise from the coupling of electrons to bosonic degrees of freedom, such as phonons. Recent 
experiments have found evidence for kinks in Ni(llO) 147 , which may be due to the electronic 
mechanism presented here. 
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6 Summary 

The combination of dynamical mean-field theory (DMFT) with LDA/GGA electronic structure 

calculations has led to important advances in our imdcrstanding and in the quantitative de- 
scription of correlated electron materials. For example, by including the p-d hybridization into 
the LDA+DMFT scheme it is now possible to clearly distinguish between Mott-Hubbard and 
charge-transfer insulators. This allowed for the computation of the local correlations and hole 
doping in NiO. By treating the local correlations and the Ni id - 02 p hybridization on the 
same footing we were able to provide a description of the full valence-band and conduction-band 
spectra of charge-transfer systems with strong hybridization. These results provide a quantita- 
tive explanation of the experimental photoemission, inverse photoemission, and ARPES data. 
Our results show that, provided the p-d hybridization is explicitly included, the LDA-I-DMFT 
is able to treat the late transition-metal oxides and charge-transfer systems in general. 

In a second application we investigated the metallization and magnetic-moment collapse 
in MnO and Fe203. Our results not only describe the experimentally observed simultaneous 
collapse of the magnetic moment and the volume in MnO but also explain the origin of the 
collapse of the magnetic monicint which originates from the increase of the crystal field splitting 
and not from the change of the band width. Since similar results were obtained for hematite 
(Fe203) we may conclude that the transition in this compound at about 50 GPa may be 
described by an electronically driven volume collapse. 

Our formulation of GGA-I-DMFT in terms of plane- wave pseudopotentials provides a robust 
computational scheme for the investigation of complex correlated materials. We applied this 
novel approach to explain various electronic and structural properties of the two prototypical 
Jahn- Teller materials KCuFa and LaMnOs in their correlated paramagnetic phase. Our results 
for the equilibrium cooperative JT distortion, the lattice constants, and the antiferro-orbital 
order agree well with experiment. This computational scheme will make it possible to investi- 
gate correlation induced structural transformations, shifts of equilibrium atomic positions, and 
changes in the lattice structure, and to perform the structural optimization of correlated para- 
magnetic solids. This opens the way for fully microscopic investigations of lattice instabilities 
which are known to occur in the vicinity of correlation-induced metal-insulator transitions. 

Finally, we offered a comprehensive explanation of kinks in the effective dispersion of corre- 
lated electrons in correlated materials with a pronounced three-peak structure. In particular we 
found that Landau quasiparticles exist only in a remarkably narrow region around the Fermi 
energy. For higher energies, but still in the central maximum of the spectral function, the ex- 
citations are described by a different dispersion relation which may be calculated analytically. 
The transition from one to the other regime is relatively sharp and leads to the kinks in the 
dispersion relation. One of the most astonishing results of this investigation is the fact that it 
is possible to describe properties of the system outside the Landau-Fermi-liquid regime fully 
analytically. 

In view of the rapid developments of the LDA-I-DMFT scheme during the last couple of 
years it is clear that this approach has a great potential for further developments. It is al- 
ready foreseeable that the LDA-I-DMFT approach will even be able to explain and predict the 
properties of complex correlated electron materials. 
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